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ABSTRACT 

This  paper  describes  a  finite  element  approach  to  quasi-TEM  analysis  of  several  different 
types  of  isolated  and  coupled  microwave  transmission  lines.  Both  the  first  and  higher-order 
ordinary  elements,  as  well  as  singular  and  infinite  elements  are  used  to  solve  for  the  potential  and 
field  distributions  in  the  cross-section  of  the  line.  Next,  the  cross-sectional  field  distribution  is 
inserted  in  a  variational  expression  to  compute  the  capacitance  per  unit  length  of  the  line  and  the 
effective  permittivity  and  characteristic  impedance  of  the  line  are  obtained  from  the  capacitance 
value.  A  perlurbational  approach  is  developed  for  estimating  the  losses  due  to  conductor  and 
dielectric  dissipation  and  computing  the  attenuation  constant. 

Lines  treatable  by  this  method  may  contain  an  arbitrary  number  of  arbitrarily  shaped 
conductors,  including  a  system  of  conductors  either  placed  above  a  single  ground  plane  or  between 
two  parallel  ground  planes,  and  inhomogeneous  dielectric  regions  that  can  be  approximated  locally 
by  a  number  of  homogeneous  subregions. 

The  results  obtained  using  the  finite  element  procedure  have  been  compared  for  various  types 
of  microwave  transmission  lines  and  have  been  found  to  agree  well  with  available  theoretical  and 
measured  data. 
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I.  INTRODUCTION 


The  accurate  prediction  of  the  characteristic  impedance,  attenuation,  coupling,  crosstalk,  etc., 
in  microstrips,  striplines  and  similar  transmission  lines  is  important  in  microwave  and  millimeter 
wave  integrated  circuits,  digital  circuit  design,  communication  and  other  applications.  The  objective 
of  this  paper  is  to  consider  microwave  transmission  lines  with  rather  arbitrary  configuration  and  to 
present  a  computer-aided  analysis  that  allows  simple  and  accurate  calculation  of  its  design 
parameters. 

In  the  past,  microwave  transmission  lines  have  been  thoroughly  investigated  by  many  authors 
and  many  approaches  to  analyzing  them  have  been  devised,  e.g..  the  Green  s  function  techniques 
[l]  -  [6].  conformal  mapping  [7]  -  [9].  variational  methods  [lO).  [ll],  Fourier  transform  method 
[12].  [l3].  Fourier  integral  method  [14],  spectral-domain  method  [l5]  -  [l7].  boundary  element 
method  [iS].  (l9l  and  finite  element  method  [20],  .All  of  the  above  methods,  with  the  e.xception  of 


the  last  two  are  restricted  in  their  application  to  transmission  lines  in  which  the  conductors  are  thin 


strips  or  the  dielectrics  inserts  base  planar  interfaces  or  both  restrictions  apply.  In  contrast,  the 


finite  element  method  (FEM)  is  capable  of  handling  transmission  lines  with  rather  arbitrary 


configurations,  since  the  lines  treatable  by  this  method  may  contain  an  arbitrary  number  of 


conductors  of  arbitrary  shape  and  inhomogeneous  dielectric  regions  that  can  be  approximated 


locally  by  a  number  of  homogeneous  subregions.  Because  of  the  generality  of  the  FF.M  approach 


;%yv-2 


the  finite  element  method  is  emploced  in  this  paper  for  the  anaKsis  of  several  representative 


microwace  transmission  lines  of  practical  interest.  The  anaKsis  is  based  on  a  quasi-TEM  model 


.cliich  IS  olten  adequate  lor  mKrowave  treqiiencies  in  communication  applications  and  for  typical 


f. 

C,*.-  ’ 


pulse  rise  limes  oi  interest  in  higii  -peed  digital  circuit  design. 


r 


H.uh  the  lirsi-  and  higher->  rder  ordinarc  elements,  as  well  as  singular  and  in!'i.''iie  elements. 


are  used  in  tlie  FL'M  algorithm  to  si)Ke  for  the  quasi-static  potential  and  the  corresponding  field 


distribution  in  a  microwace  transmission  line.  The  capacitance  per  unit  length  of  the  line  is 


obtained  from  a  variational  e\pres-uin.  anil  the  elTeclive  dielectric  permittivity,  characteristic 
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impedance  and  phase  velocity  are  also  calculated.  Finally,  a  perturbation  method  is  used  to 
compute  the  losses  due  to  both  the  conductor  and  dielectric  dissipations. 

The  characteristic  impedance  and  loss  characteristics  have  been  calculated  for  various 
transmission  lines  of  interest  and  very  good  agreement  with  available  theoretical  and  experimental 
data  has  been  obtained. 

n.  FINITE  ELEMENT  MODEL  OF  A  MICROWAVE  TRANSMISSION  LINE 

Consider  a  microwave  transmission  line  with  an  arbitrary  cross  section  consisting  of  a  number 
of  arbitrarily  shaped  conductors  and  inhomogeneous  dielectric  regions  which  can  be  approximated 
locally  by  homogeneous  subregions  (Fig.  1).  Let  us  assume  that  the  line  is  uniform  along  its 
longitudinal  (2)  axis  and  let  €  and  /j.  (=  /j.o)  denote  the  permittivity  and  permeability,  respectively, 
of  the  medium  in  each  homogeneous  subregion. 

We  assume  a  quasi-TEM  model,  i.e..  that  the  dominant  mode  propagating  along  the  line  is  a 
TE.M-mode.  Under  this  approximation,  the  problem  is  reduced  to  that  of  finding  the  scalar 
potential  function.  0.  .satisfying  the  Poisson  s  equation 

A0  =  -^.  (1) 

€ 

.subject  to  the  (appropriate)  boundary  conditions 

0  =  0,,  on  r  1  =  0  on  r^.  (la) 

where  p  is  the  volume  charge  distribution,  e  is  the  dielectric  permittivity  and  Fi  and  r2  are 
complementary  portions  of  the  boundary  f  of  the  domain  S  in  which  Eqn.  (1)  is  defined. 

In  order  to  apply  the  finite  element  approach  to  the  problem  at  hand  we  formulate  it  in 
\  ariational  terms.  The  correct  solution  of  I’qn.  ( I )  is  one  that  minimizes  the  energy  functional 

F  =  ^  f  6  V0  i  -  dS  -  f  (xAdS  (2) 

Eqn.  (1)  IS  the  Euler  equation  of  the  functional  E  (27)  satisfying  the  es.seniial  boundary  condition 
of  Dirchlei  type,  v  12.. 


and  the  natural  boundary  condition  of  Newman  type.  i.e.. 

=  0.  on  r-).  (2b) 

dn 

■According  to  the  FEM  procedure,  the  cross  section  of  the  line,  domain  S.  is  subdivided  into 
finite  elements  in  an  arbitrary  manner  provided  that  all  the  dielectric  interfaces  coincide  with  the 
element  sides.  .Although  a  variety  of  different  elements  can  be  chosen,  the  triangular  first-  or 
higher-order  [21]  elements  are  adopted  in  this  study-  It  has  been  shown  [2l]  that  the  accuracy  and 
efficiency  of  computation  can  be  substantially  enhanced  via  the  use  of  the  high-order  elements  as 
compared  to  the  case  where  only  the  first-order  elements  are  employed,  if  the  same  number  of 
nodes  is  used.  On  the  other  hand,  a  large  number  of  simple  elements  is  clearly  advantageous  where 
a  complicated  boundary  shape  needs  to  be  modeled. 

However,  when  there  are  field  singularities  caused  by  edges  in  the  cross  section  of  the  line,  a 
very  fine  mesh  of  first-  or  high-order  elements  is  required  to  obtain  an  accurate  solution.  However, 
in  order  to  improve  accuracy  and  reduce  the  number  of  necessary  nodes,  singular  elements  [22]  are 
u.sed  in  this  study.  .A  singular  element  with  labeled  nodes  and  associated  triangular  polar 
coordinate  system  (p.  <t)  is  shown  in  Figure  2. 

The  global  Cartesian  coordinates  a.'e  then 

X  =  .t  1  +  p  [(.t;;  —  .t  l)  +  cr  (.V  1  —  V;)]  Ou) 

y  =  y  ,  +  p  [(y  :  -  y  I )  +  O-  (y -  y  , )] 

The  scalar  potential  distribution  function  that  takes  into  account  the  field  singularity  at  the 
node  1  (p  =  0)  has  the  following  form: 

6  =  (t*!  i  1  —  p' '  -I-  <I>;p^  (I  —  cr  ■  'P^p'cr.  i.  <  \  <  1.  (3b) 

where  coefficient  \  is  chosen  in  accordance  with  the  Meisner  edge  condition  [23].  and  <l>i.  4>2-  *^3  ^re 
corresponding  nodal  potentials.  The  linear  cariation  of  the  scalar  potential  along  the  side  2-3  is 
provided  so  that  the  element  is  compatible  w  ith  first-order  ordinary  elements. 


Although  open-type  transmission  lines  can  be  treated  by  the  conventional  finite  element 
method  with  the  shielding  far  away  from  the  region  of  interest,  the  computation  efficiency  can  be 
substantially  improved  if  infinite  elements  [24]  are  used.  Consider  first  the  line  consisting  of  a 
number  of  conductors  arbitrarily  placed  between  two  parallel  ground  planes  (Fig.  3a). 

The  entire  domain  can  be  divided  into  the  near-field  region  (n.f.).  which  is  the  region  of 
interest  and  the  far-field  region  (f.f.)  which  is  unbounded.  They  have  a  common  boundary 
referred  to  as  the  far-field  boundary  (f.f.b).  The  near  field  (n.f.)  region  is  then  divided  into  finite 
triangular  elements  in  the  usual  manner,  while  the  f.f.  region  is  divided  into  infinite  elemenLs. 
Each  infinite  element  has  two  common  nodes  (lying  on  f.f.b)  with  an  ordinary  first-order  element 
and  two  sides  parallel  to  the  x-a.xis. 

Consider  the  infinite  element  with  nodes  1  and  2  (Fig.  3b)  and  introduce  normalized 
coordinates 


I  =  Ti  =  __L 

'i  .''.--vi 


.  .t  ,  <  <  CO.  V  I  <  V  <  V  -. 


(4a) 


Since  in  the  first  approximation  the  far-field  of  the  system  is  equivalent  to  a  dipole  field,  the 
electric  scalar  potential  within  the  element  is  chosen  in  the  following  way: 


d,  r.  [  (I)l(  l-T))  + 


(4b) 


here  and  are  potentials  of  nodes  1  and  2.  respectively. 

The  next  tv  pe  of  open  region  problems  is  a  system  of  conductors  above  a  ground  plane.  The 
procedure  is  similar  to  that  m  the  previous  ca.se.  The  entire  domain  is  divided  into  n.f.  and  f.f. 
regitins  (Fig.  4a)  The  n.l.  .-egion  is  die  ided  into  the  usual  triangular  mesh  and  the  f.f.  region  is 
dr.  i.led  into  infinite  ele.ments  v'f  two  different  types.  I  and  II.  The  infinite  element  l\  pe  1  is  the 
element  preciously  discussed.  C>)nsider  now  the  infinite  element  of  type  II  with  nodes  1.  2  and 
radial  sides  intersecting  at  point  (  v  y  }  i  I'lg  4b)  and  introduce  triangular  polar  coordinates  p.c 


cc  Inch  are  related  to  the  global  ('artesian  coordinates  bv  the  relations 


GROUND  PLANE 


r  igure  3a.  !  me  con.-^isun^  ot  a  sysiem  of  conductors  placed 
bet'A,een  fa.o  ground  planes. 
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X  =  t,  +  p  [  (.r  ,-t,  ,  +  cr(jt,-x  ,)  ] 
y  =  y  +  p  [  (  y  y,,  )  +  (T(y  y  , )  ] 


The  scalar  poieniial  disinbution  within  the  element  is  approximated  by 

d>  =  —  [  0,(1— cr)  +  OiO-  ] 

P 


where  O,.  Oi  are  nodal  potentials,  for  element  type  lla  v„  .  y,  =0 

For  both  types  of  infinite  elements,  the  linear  variation  of  the  potential  function  along  the 
side  1-2  (common  side  w'ith  a  corresponding  ordinary  element)  is  provided.  These  elements  are 
compatible  with  ordinary  first-order  elements  and  satisfy  completeness,  finiteness  and  radiation 
conditions  [24], 

Once  the  potential  distribution  is  known,  the  capacitance  per  unit  length,  (b.  is  readily 
obtained  from  the  cariaiional  expression. 


VC  here  W  is  the  energy  per  unit  length  of  the  line  and  V  is  the  conductor  potential.  Then  the 
effective  (dielectric)  permittiv  ity  is 


F*  •  *  -  » 


c; 


where  C  is  the  capacitance  of  the  actual  line.  C,  is  the  free-space  capacitance  of  the  line.  The  phase 
velocitv  IS 


lere  c  i>  il-,e  .eli'c:'  .  .1  liclil  in  I  ,ve  '•p.i».e  I  he  v. Iiaraclerisi  k  impedance  is 


,> 
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The  expressions  given  in  (6a)-(d)  are  for  a  single  line,  but  the  same  procedure  can  be  applied 
to  find  the  capacitance  and  characteristic  impedance  for  even-  and  odd-modes  on  two  coupled  lines 
or  to  find  capacitance  and  impedance  matrices  for  a  system  of  n-lines.  The  electromagnetic  field 
distribution  can  also  be  obtained: 


=  fi.  =-  -E  ^2 

V 

VC  here  i  is  unit  vector  in  the  z  direction  and  V  is  the  intrinsic  impedance  of  the  medium 


(7a) 


V  -  (Mo  )'■  (7b) 

Next,  we  employ  a  perturbational  approach  to  solve  for  the  attenuation  constants  due  to  dielectric 
and  conductor  losses 


a, 


p.,  r. 

2p  ■  ~  2P 


(Sa) 


where  P.  is  ilie  ume-averaged  power  tlow  along  the  line,  and  Pj  and  P,  are  the  lime-averaged 
povcers  dissipated  in  the  dielectrics  .ind  conductors,  respectively. 

Dielectric  los.ses  are  calculated  using  the  formula 
P .  =  weianfi  J  E  -JS  (Sb) 


VC  here  the  loss  tangeni.  lan  5,  is  .issumed  to  be  sufficiently  small  so  that  the  perturbed  fields  can  be 
.ipproximated  b\  the  tields  for  the  lossless  condition  E.  . .  w  =  2iZI  is  the  angular  frequency 
.ind  ,S  IS  the  area  ol  ^ross  section  covered  by  the  dielectric. 

I  '-ses  due  to  •he  imperic.l  ^>'ndiui.-!N  are  obtained  v  i.i  the  v-  nventional  perturbation 

;  I'ri'.iul.i  i2-' 

/’  =  R  J  II  y,  dl 

.c  here  A’  is  the  surface  resistance,  .uid  II  ,  .  is  the  magnitude  of  the  tangential  magnetic  field 
a’  ih.e  ^.'iiducting  surlaces  t  line  ('  ivU'  the  lossless  case. 


.  .S 


The  average  power  propagating  along  the  line  is  given  by 

=Ref  (E,  dS  (8d) 

A' 

where  S  js  the  complete  cross  section  of  the  line. 

For  the  conductor  loss  calculation,  special  finite  element  matrices  have  been  derived  [26]. 
Consider  a  high-order  element  where  vertices  are  labeled  as  nodes  i.  j.  k  (Fig.  5). 

Suppose  that  the  side  j-k  is  an  imperfect  conductor  of  surface  resistance  ■  Then  the  losses  in  that 
side  can  be  expressed  as 

K 

P'.  =  — —  P,  coc  B,  +  Pj  cot  Bj  +  Pi^  cot  4> 

T 

where  [O]  is  a  column  vector  of  nodal  potentials  for  the  element,  B, .  B^  .  are  the  included  angles 
at  vertices  i.  j.  k.  respectively  and  P, .  Pj  .  Pi  are  matrices  of  the  numerical  coefficients  which  do 
not  depend  on  the  triangle  shape. 

in.  NUMERICAL  AND  EXPERIMENTAL  RESULTS 

On  the  basis  of  described  finite  element  procedures,  two  different  computer  software  packages 
have  been  developed  for  solving  the  problem  at  hand.  The  first  of  these  is  employed  for 
semiautomatic  mesh  generation  of  first-  or  high-order  ordinary  elements  including  singular  and 
infinite  elements  if  necessary.  The  second  package  calculates  the  potential  and  field  distributions, 
capacitance  per  unit  length,  characteristic  impedance,  effective  permittivity,  and  attenuation  due  to 
conductor  and  dielectric  losses.  The  linear  system  of  equations  derived  from  the  application  of  the 
FEM  method  is  rather  spar.sc  and  special  techniques  for  the  solution  of  such  equations  can  be 
employed  to  achieve  enhancei!  computational  efficiency.  Two  such  techniques  are  direct  envelope 
and  band  method,  with  node  reordering,  which  are  u.seful  for  smaller  problems  with  a  few 
hundred  nodes.  For  large  systems  of  equations  (a  few  thousand  nodes),  the  conjugate  gradient 
method  i.s  suggested  because  it  requires  the  storage  of  only  non-zero  elements.  Preconditioning  is 


accomplished  via  incomplete  Cholesky  decomposition  [27]  which  leads  to  the  acceleration  of  the 
Iteration  procedure. 

To  demonstrate  the  quality  and  the  accuracy  of  the  finite  element  method  analysis  of  the 
previous  section,  the  solutions  for  sample  problems  are  given  and  are  compared  with  available 
theoretical  and  experimental  data,  as  well  as  with  measured  results. 

The  first  example  considered  is  the  problem  of  a  shielded  multilayered  thick  strip 
transmission  line  shown  in  Fig.  6.  The  corresponding  triangular  mesh  for  one  half  of  the  line  is 
presented  in  Fig.  7.  The  characteristic  impedance  of  the  line  is  calculated  for  a  few  different 
thicknesses  of  the  strip,  t.  b.  as  a  function  of  ratio  2a/b  (width  and  height  of  the  shielding)  and  is 
represented  by  the  .solid  line  m  Fig.  h.  Corresponding  results  obtained  by  the  Green’s  function 
method  [3]  are  denoted  by  the  dashed  line  and  good  agreement  between  the  two  sets  of  results  is 
evident.  It  is  found  that  the  influence  of  the  side  walls  on  the  characteristic  impedance  is  negligible 
when  the  walls  are  sufficiently  removed  from  the  strip.  2a''b>5. 

The  next  example  is  an  inverted  microstrip  whose  generic  cross  section  (see  [17.  19])  is  shown 
in  Fig.  9.  Calculated  results  for  the  characteristic  impedance  as  a  function  of  aspect  ratio  2w/h. 
denoted  by  the  solid  line,  are  compared  wi'h  the  numerical  results  obtained  via  the  boundary 
element  [BEM]  method  [l9],  shown  by  the  dashed  line  m  Fig.  10  Good  agreement,  within  5%.  is 
obtained.  Data  tor  v of  the  line  as  a  function  of  aspect  ratio  2w  'h  are  presented  in  Fig.  11. 
1  mite  element  results  are  denoted  b>  solid  lines  and  B[:M  results  by  dashed  lines.  The  measured 
points  shoNcn  m  the  same  figure  \cere  obtained  via  time-domain  reflectometer  measurements  [28]. 
It  can  be  noted  that  FEM  results  are  «. loser  to  the  experimental  results  than  the  BEM  data.  For  the 
same  example,  attenuation  constants  itue  to  conductor,  q  .  and  dielectric.  Qj  .  lo.sses  are  calculated 
as  (unctions  of  aspect  ratio  2VC  h  ami  shoccn  in  Fig.  12  as  solid  curves.  These  data  are  compared 
cciih  BEM  [id]  and  spectral-domain  method  (l7]  results  shown  in  the  .same  figure  and  a  very  good 
agreement  is  found. 
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li  should  be  pointed  out  that  the  accuracy  of  the  calculation,  especially  of  the  loss 
calculations,  can  be  significantly  enhanced  via  the  use  of  singular  elements.  The  odd-  and  even- 
modes  propagating  on  coupled  lines  can  be  treated  by  the  same  finite  element  procedure.  Some  of 
the  results  for  coupled  lines  are  presented  in  [29]. 

In  order  to  verify  the  accuracy  of  the  finite  element  procedures,  the  characteristic  impedance 
and  losses  of  a  shielded  microsirip  line  have  been  measured.  The  experiment  has  been  done  on  the 
line  whose  cross  section  and  the  dimensions  are  shown  in  Fig.  13.  By  using  a  HP  8510  Network 
.Analyzer,  time-domain  reflectometer  mea.surements  of  characteristic  impedance  (Fig.  14).  and 
frequency-domain  measurements  of  the  Sn  parameter  of  open-  and  short-circuited  lines  (Fig.  15) 
have  been  done.  The  measured  value  of  characteristic  impedance  is  48  ohms  while  the  calculated 
value  is  48.1  ohms.  Excellent  agreement  between  measured  and  calculated  values  of  the  Sn 
parameter  of  the  line  is  found. 

IV.  CONCXUSIONS 

In  this  study,  a  finite-element  procedure  for  quasi-TEM  analysis  of  microwave  transmission 
lines  is  presented.  The  method  is  capable  of  handling  microwave  transmission  lines  with  rather 
arbitrary  configurations.  Both  first-  and  high-order  triangular  elements  are  employed.  In  the  case 
of  open  region  problems  (conductors  over  one  or  between  two  parallel  ground  planes)  and  edge 
singularity  problems,  accuracy  and  efficiency  of  the  method  are  enhanced  via  the  use  of  infinite  and 
singular  elements.  For  the  solution  of  finite  element  linear  equations,  the  band  or  envelope  matrix 
method  (with  reverse  ('uthill-McKee  reordering)  is  used  for  smaller  systems  and  the  incomplete 
Cholesky-conjugate  gradient  meihiid  for  large  systems. 

Computer  software  for  semiautomatic  mesh  generation  as  well  as  for  calculation  of 
capaciL  ince  per  unit  length,  elfectice  permittivity,  characteristic  impedance  and  attenuation  due  to 
conductor  and  dielectric  losses  has  been  developed. 

Results  obtained  for  various  tv  pes  of  microwave  transmission  lines  are  in  very  good 
agreement  with  available  theoretical  and  e.xperimental  data. 


Sii  Re 

Sll 

1 og  MAG 

REF  0.0  Units 

REF 

0.0  dB 

i  5.0  mUn i t  s/ 

A 

10.0  dB/ 

V  4.6911  mU. 

4 

29.749  dB 

50  GHM  LINE  MAY  21 

fma 

111 


5uir  “ 


START 

STOP 


150 .0  p  s 
1.2  ns 


START  0.045000000  GHz 
STOP  18.045000000  GHz 


f  igure  14.  Time-domain  measured  results. 
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Figure  15.  Frequency -domain  measured  results  of  parameter 
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